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Phase diagram of magnetic domain walls in spin valve nano-stripes 
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We investigate numerically the transverse versus vortex phase diagram of head-to-head domain walls in 
Co/Cu/Py spin valve nano-stripes (Py: Permalloy), in which the Co layer is mostly single domain while the 
Py layer hosts the domain wall. The range of stability of the transverse wall is shifted towards larger thickness 
compared to single Py layers, due to a magnetostatic screening effect between the two layers. An approached 
analytical scaling law is derived, which reproduces faithfully the phase diagram. 



The fundamental study of magnetic domain 
walls (DWs) has found a playground in lithographically- 
defined nano-stripes. In the range of widths from 
a few tens to a few hundreds of nanometers their 
complexity is dramatically reduced with respect to 
extended thin films, while still retaining a few internal 
degrees of freedom leaving a rich physics. For stripes 
with in-plane magnetization, head-to-head DWs are of 
transverse (TW) or vortex (VW) type (FIG. Ib)^'^, and 
may be characterized by a chirality, asymmetry and/or 
polarity. Their propagation under the stimulus of a 
magnetic field or a current of spin-polarized electrons 
is intrinsically precessional and reveals effects such 
as DW inertia^ and a so-called Walker limit beyond 
which periodic DW transformations occur '^^. Due to 
their small size and fast dynamics, magnetic domain 
walls (DWs) are promising candidates in the area of 
information processing and storage^'^. 

Most reports so far have considered DWs in a sin- 
gle layer of magnetically-soft materials, such as Permal- 
loy (Py) or CoFeB alloys. Magnetic trilayers F1/NM/F2 
(with Fi a ferromagnet, and NM a non-magnetic ma- 
terial) also deserve attention as they play a key role in 
spintronics, especially in devices using Giant Magneto- 
Resistance. These stacking are called spin valve for a 
metal spacer layer, and pseudo spin valve or magnetic 
tunnel junction for an insulator spacer layer. Fl is a soft 
layer in which the DW moves, while F2 is a reference 
layer that is supposed to remain uniformly magnetized. 
Unusually-large current-driven DW mobilities with speed 
exceeding 600 m/s have been reported in such stacks^'^. 
Tentative explanations point at either out-of-plane spin 
accumulation and torque^^'^^. Oersted field^^ or magne- 
tostatic coupling between the layers that may modify the 
statics or dynamics of the DW^^'^^. 

While the role of magnetostatic coupling is well estab- 
lished in thin films with the formation of quasi-walls in 
the supposedly uniformly-magnetized layer ^^, it has been 
scarcely addressed ^'^ in stripes. In this Letter we deter- 
mine the phase diagram of head-to-head DWs in spin 
valves with in-plane magnetization (TW versus VW), 
combining numerical simulation and analytical modeling. 
For this purpose, we consider Co/Cu/Py nano-stripes 
with various widths w and thickness of the individual lay- 



ers tco, ^Cu and tpy. The DW is located in the Py layer, 
while the Co layer is initialized with uniform magnetiza- 
tion. We predict that the stability of TWs is enhanced 
towards larger Py thickness compared to single layers^ '^, 
due to the screening in the Co layer of the stray field 
arising from the DW in Py. 

We used finite- differences micromagnetic codes, both 
OOMMF^^ and our home-made GL-FFT code^^. The 
cell size in all simulations presented here is 4 x 4 x 5 nm^ ; 
fo significant difference is found for a cell size 4 x 4 x 
2.5 nm^. Effects of the finite length of the stripe in the 
simulations were reduced by either fixing at both ends the 
magnetic moments parallel to the stripe direction, or by 
compensating the surface magnetic charges at both ends, 
however with no constraint on magnetization. The length 
L of the stripes was set equal or larger than about 20w, 
for which finite-size effects are identically small under ei- 
ther of the above procedures. Material parameters are 
MoMco = 1.7593 T, Aco = 30pJ/m, fioMpy = 1.0053 T, 
Apy = lOpJ/m, and zero magnetocrystalline anisotropy 
in both layers. No exchange coupling between Py and 
Co is considered either. Different initial conditions were 
chosen in Py to set a TW or VW in the Py layer, while 
magnetization in Co is set uniform along the stripe di- 
rection; both layers are then set free to evolve during 
the simulation. We considered head-to-head DWs (two 
domains facing each-other, see FIG. la) with no loss of 
generality, as tail-to-tail DWs are equivalent under time- 
reversal symmetry (M -^ — M). Like for single layers, 
TWs and VWs coexist as (meta) stable states in a large 
range of geometrical parameters. The iso-energy curve 
separating the areas where either a TW or a VW is low- 
est in energy, was determined for each fixed w through a 
parabolic interpolation of the energy difference for three 
different values of tpy. We checked that this procedure 
yields a phase diagram of the single Py layer in very close 
agreement with existing reports^'-^. 

We first discuss qualitatively the difference between a 
DW in a single layer and in a spin valve. The head- 
to-head DW in Py is associated with a total magnetic 
charge -\-2tpywMs (FIG. la). While VWs are similar 
in both systems, clear differences are found for TWs. 
In a single layer, TWs are symmetric for small thick- 
ness and width (FIG. 2). Going towards a larger thick- 
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FIG. 1. (a) Cross-sectional sketch of a head-to- head domain 
wall and the associated charge in a spin valve. Plane views 
of a vortex wall and transverse wall in a 300 nm-wide (b) sin- 
gle layer of Py(5nm) and (c) a Co(5 nni)/Cu(5 nm)/Py(5 
nm) spin valve, (d) thickness-integrated maps of magnetiza- 
tion from (c), see Eq. (1). The red (blue) color indicates the 
positive (negative) transverse component of magnetization, 
normalized to the magnetization of the layer under focus in 
(b) and (c), and to half the magnetization of Co in (d). 



ness or width, the TW becomes asymmetric, reduc- 
ing the magnetostatic energy by spreading the magnetic 
charges^ (FIG. lb); this corresponds to the onset of zig- 
zag walls in continuous films. The transition from sym- 
metric to asymmetric TW is of second order (contin- 
uous), so that for a given geometry only one type of 
TW exists, either symmetric or asymmetric. Asymmetric 
TWs are found for a large range of values of width and 
thickness, although they are the state of lowest energy 
(over the VW) only in a narrow range of values (FIG. 2). 
In a spin valve the Co layer deviates from uniform mag- 
netization due to the stray field arising from the DW 
in Py, which in turn creates a stray field acting on the 
Py layer. This magnetic screening effect lowers the en- 
ergy compared to a single layer, as already pointed out^^. 
The stray field arising from the DW in Py is parallel 
or antiparallel to the initial direction of magnetization 
of Co, depending whether the left or right side of the 
DW (FIG. la) is considered. Thus, the magnetic config- 
uration of Co is expected to be asymmetric, and so will 
be the stray field arising from Co and acting on the Py 
layer. A TW in a spin valve is therefore expected to be 
always asymmetric by nature, due to the unidirectional 
magnetization in the Co layer. 

FIG. Ic shows a VW and a TW in a spin valve, with the 
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FIG. 2. Phase diagram of head-to-head domain walls in 
stripes, with the boundary between symmetric and asymmet- 
ric TWs (A), and TW and VW (o) in a single layer, and 
asymmetric TW and VW in a spin valve (•). From bottom 
to top, the three lines are: a 1/tpy fit to the TW-VW data for 
a single layer, a thickness translation of the former to fit the 
TW-VW data for a spin valve, and another thickness trans- 
lation with tsh{tcu = 0), see Eq. (3) and text for details. 



same stripe width and Py thickness as in FIG. lb. The 
magnetic screening in Co is clear from the non-uniformity 
of magnetization. As expected from the symmetry argu- 
ments given above the TW is asymmetric, to a larger 
extent than in the single layer case^^. Both the VW and 
TW have a larger width, associated with the emergence 
of tails. Indeed, as a DW profile results from the balance 
of magnetostatic with exchange energy, the decrease of 
magnetostatic energy allows to reduce exchange energy 
via the increase of the DW width. This is consistent with 
simple analytical models ^^ derived for charged (head-to- 
head) domain walls in extended spin- valve thin films. 

The flux closure taking place between Py and Co can 
also be illustrated by the map of magnetization Mint inte- 
grated over the two layers along their normal (FIG. Id): 



Mint = 



fM{z)dz (1) 



This map shows similarities with a VW, which highlights 
the principle of reduction of the total energy: the flux is 
better closed than in a single layer TW, while avoiding 
the cost associated with the vortex core in a single layer 
VW. In this view the DW tails of the individual layers 
are also absent. 

As the share of magnetostatic energy is larger in a TW 
than in a VW^, the inter-layer closure of flux is lowers 
the energy of the TW more than that of the VW. This 
suggests that the TW should be the ground state in a 
larger range of geometrical parameters in a spin valve 
than in a single layer. This is confirmed by a set of sim- 
ulations with varying tpy and w, illustrated in FIG. 2 for 
tco = ^Cu = 5nm. Notice that in a single layer TWs 



with opposite asymmetries (left or right) are degenerate 
in energy, while this degeneracy is lifted in a spin valve, 
due to the unidirectional magnetization in the underlying 
Co. In the phase diagram (FIG. 2) we considered only 
the TW with the lowest energy. 

We now derive a simple analytical model to grasp the 
main features of this diagram. McMichael and Donahue 
already noticed that the iso-energy line in a s ingle la yer 
follows the power law tw = C A^ , with A^ = ^J AjK^ the 
dipolar exchange length, K^ = {l/2)fioM^ the dipolar 
constant, and C a constant. This law can be derived 
qualitatively by balancing the energies at play in each 
domain wall^. These authors argued that the energy of 
the TW can be estimated from the lateral demagnetizing 
coefficient scaling like t/w^ and the volume of the TW 
of the order of tw^ (thus t^w as a whole), while that 
of the VW is related to the energy of the vortex core, 
scaling with tA'^K^. The exchange energy in the ^ 90° 
Neel sub- walls can be ignored as the total length of these 
walls is identical in a TW and a VW. The numerical value 
C ^ 61 — 64 must be provided by simulations^'^. 

This scaling model can be adapted to a spin valve, 
however the above arguments must ffist be discussed in 
more detail. It shall ffist be noted that in any head- 
to-head DW, either TW or VW, most of the energy is 
of magnetostatic origin and related to the head-to-head 
charges, so that the above argument must be refined. 
Simulation and magnetic force microscopy of a head-to- 
head DW in a single layer show that the total charge Q = 
2twMs is nearly uniformly spread over the area of the 
DW^, which is ^ w^ for a TW, and ^ 2w^ for a VW. The 
associated densities of surface charges are then axw = 
2{t/w)Ms and ayw = {t/w)M^. The resulting volume 
density of stray field energy scales with <j^ and extends 
over a typical distance w above and below the stripe. 
This results in an excess of magnetostatic energy in the 
TW over the VW of ^ 2t^wK(X. We now consider a spin 
valve. In the limit tpyMpy > tco^Co? which is suitable 
for tco = 5 nm and the boundaries of the phase diagram 
of practical interest^'^, the deformation in Co reaches its 
maximum while not being enough for a full screening of 
the charges arising from the DW in Py. The above scaling 
laws can be rewritten accordingly and the following law 
is derived again with some approximations (see below), 
making also use of the partial fiux closure in spin valves 
to refine the lateral demagnetizing coefficient^^: 



w{t-Ui,)^CAl 



(2) 



This suggests that the phase diagram for a spin valve 
is based on that of a single layer, shifted towards higher 
thicknesses by a value t^h- This is in very good agreement 
with the results of numerical simulations (FIG. 2). Quan- 
titatively, the analytical formula for tgh is not simple, im- 
plying both power and logarithmic functions. However, 
thanks to the slow variation of the latter, tgh can be ap- 
proximated faithfully with tco(^Co/^Py) ^ l-75tco- 
For tco = tcu = 5 nm the analytical model thus yields 



1.75tco ^8.75 nm, while tgh determined from simulations 
is 5.4 nm. This discrepancy is linked with the magnetic 
screening in the Co layer being partial. This can be un- 
derstood as in the scaling model we neglect the magne- 
tostatic energy stored in the volume of the Cu spacer, 
locus of a magnetic field arising between the charges of 
opposite sign in Py and in Co. FIG. 3 shows t^h re- 
sulting from micromagnetic simulations as a function of 
tcu for ^ = 100 nm and tco = 5nm. For large tcu the 
problem tends towards the case of a single layer because 
the screening is no more effective. The decay is close 
to exponential, with a fit providing the extrapolation 
4h(^Cu = 0) ^ 8.60 nm. This figure is now in excellent 
agreement with the analytical model, which suggests the 
following empirical law: 



t'sh 



, Mco 
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*0 



Py 



(3) 



where to ~ 10 nm is derived from FIG. 3. Additional mi- 
cromagnetic simulations show that to increases with w, 
e.g. to = 13 nm and 16 nm for w = 300 nm and 500 nm, 
respectively. This trend can be understood as the ini- 
tial stray field (to be reduced through screening) arising 
from the DW extends outside the stripe over a volume 
scaling with w^, while the magnetic field arising between 
the charges of opposite sign in Py and in Co (a side cost 
of the screening effect) applies in a volume scaling with 
tcu'^^- Thus, the relative gain of energy through screen- 
ing decreases with decreasing the stripe width. Then, 
Eq. (3) can be applied for other Co thicknesses, or more 
generally to other couples of magnetic materials if sub- 
stituting the proper spontaneous magnetization values in 
Eq. (2), and fitting the decay of t^h like in FIG. 3. 

One must finally keep in mind that for tpyMpy < 
tcoMco a very effective screening of the head-to- head 
DW charge is possible, so that deviations from Eq. (2) 
are expected. In that case, the energy of the TW is dras- 
tically reduced. The iso-energy line becomes very fiat 
beyond this limit, an effect that starts to be visible for 
the largest widths in FIG. 2. 

To conclude, we have investigated with numerical sim- 
ulation and analytical scaling laws the phase diagram 
of head-to-head domain walls in F1/NM/F2 spin valves, 
with a domain wall in one layer and no domain wall in 
the other layer. We showed that the range of stability of 
transverse versus vortex walls is enhanced due to a mag- 
netic screening effect. The iso-energy line in spin valves is 
translated towards a larger thickness with respect to sin- 
gle layers, by a value decreasing approximately exponen- 
tially with the spacer thickness. This enhanced stability 
and larger width provide a magnetostatic contribution to 
the high mobility of DWs in trilayers, provided that the 
damping in Co is not too large^'^. Effects on the domain 
wall inertia and automotion'^ are also expected. 
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FIG. 3. Thickness shift of the iso-energy hne in a spin valve 
compared to a single layer, for it; = 100 nm and tco = 5 nm. 
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